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(N 
> 

■ Isolated isothermal spheres of N gravitationally interacting points with equal mass 

^ ' are believed to be stable when density contrasts do not exceed 709. That stability limit 

o ■ 

^ ! does, however, not take into consideration fluctuations of temperature near the onset of 
instability. These are important when N is flnite. 

O 
}^ 

' Here we correlate global mean quadratic temperature fluctuations with onset of insta- 

bility. We show that such fluctuations trigger instability when the density contrast reaches 
^ ! a value near 709 • exp(— S.SA^"^/"^). These lower values of limiting density contrasts are 
significantly smaller than 709 when N is not very big and this suggests (i) that numerical 
calculations with small N may not reflect correctly the onset of core collapse in clusters 
with big N and (ii) that a greater number of globular clusters than is normally believed may 
already be in an advanced stage of core collapse because most of observed globular clusters 
whose parameters flt quasi-isothermal conflgurations are close to marginal stability. 
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1. Introduction 

Bounded isothermal spheres of pointhke stars of equal mass interacting gravitationally 
are unrealistic but instructive models to learn about the early stages of gravitational col- 
lapse in stellar clusters. Their great advantages are mathematical simplicity and the fact 
that other factors of evolution than the "gravothermal catastrophe" like escape of stars 
to infinity, equipartition of energy among stars of different masses or formation of hard 
binaries are switched off by the model. Isothermal spheres have been studied in statistical 
mechanics (Antonov 1962, Horwitz and Katz 1978), in thermodynamics (Lynden-Bell and 
Wood 1968, Katz 1978) and in numerical experiments (Hachisu et al. 1978, Inagaki 1980). 
The various studies came up with consistent results and confirmed the following picture of 
slow evolution of isothermal spheres through different stages of global quasi-static isother- 
mal equilibrium with increasing entropy and density contrasts. It is nicely described In 
Galactic Dynamics by Binney and Tremaine (1987). 

Consider a bounded isolated isothermal sphere of radius R with N point particles 
of mass m attracting each other gravitationally. The total energy E is a function of the 
inverse global temperature T shown in figure 1 . The winding-in of pairs of thermodynamic 
equilibrium quantities like E and P = 1/kT as shown in figure 1 is also a property of 
polytropic gas spheres (Chandrasekhar 1934), cold white dwarfs (Harrison et al. 1965), hot 
isentropic stars (Thorne 1966) and relativistic gases with P/p =constant (Chandrasekhar 
1972). Each point on the E{P) line, like A, B, C, etc. is associated with an equilibrium 
configuration. The branch of thermodynamically stable equilibria indicated on figure 1 
ends at point D where the density contrast is ~ 709 and the energy is also the smallest 
for equilibrium configurations (Antonov 1962; Lynden-Bell and Wood 1968). All points of 
density contrasts bigger than ~ 709 represent unstable configurations (Katz 1978). 
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Evolution is down the slope in figure 1 say from point A towards point D, towards 
higher entropy and higher density contrasts. What happens to a slowly evolving isothermal 
sphere when it becomes unstable? Various interesting thought experiments have been 
performed and vividly described in Lynden-Bell and Wood and in Binney and Tremaine. 
The system becomes unstable when it reaches point D and a "gravothermal catastrophe" 
develops. Notice that from B to D the entropy is a local maximum, while the heat capacity 
Cv < since the slope is positive. From D to F the entropy extremum is a minimum but 
Cv > 0. 

Notice also that those thought experiments are based on the assumption that the 
system reaches a new local entropy minimum in a finite time. In fact the system spends a 
limited time at each point of the linear series in (quasi-) equilibrium. It is thus assumed that 
during the evolution, fluctuations, and in particular mean global temperature fluctuations 
(which we shall precisely define in a moment) are small, i.e. that the probability of a 
big fluctuation is negligibly small. This assumption is certainly not correct near point D 
where fluctuations of global temperature become infinite (Horwitz and Katz 1977 referred 
to below as HK77). One of the main objects of the present work is to try to assess how 
close to point D the evolution can proceed before the system becomes unstable as a result 
of big fluctuations. If the system spends an unlimited time in a stable equilibrium with 
energy less negative than Ed-, it will eventually evolve to a non isothermal state because 
big fluctuations of temperature have a chance to develop and bring the system to a state 
of entropy which would increased indeflnitely (see below) . 

Now a few words about temperature fluctuations in isothermal spheres which are at the 
center of our considerations. The global temperature T is a global equilibrium parameter 
but it has also a local meaning (see Lynden-Bell and Wood again) ; T is related to the mean 
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square velocity by the same relation as in an ideal gas \kT = ^mv^. The mean square 
velocity is that of the stars which happen to be near a point at a given moment. At later 
times stars left the vicinity and other stars came around. A local temperature can also 
be defined if the cluster is not in global thermal equilibrium but such a local temperature 
would differ from point to point and from time to time. A global equilibium temperature 
does then not exist. The proper method to tackle the problem is then in kinetic theory 
not statistical mechanics. 

However, if an out of equilibrium entropy can be defined in statistical physics for a 
system with a given energy E [see HK 1977], say w, than a global out of equilibrium 
temperature T can be defined as well by dw/dE = 1/kT = p. That T will clearly be 
equal to T at the point where w is equal to its maximum S and fluctuations of T can 
be calculated like in classical statistical physics (see Landau and Lifshitz 1985). Thus we 
may state without proof (but we shall give the prove below) that the mean quadratic 
fluctuations of the global temperature (AT)2/T^ (from now on temperature will always 
mean global temperature) is, like in classical thermodynamics, of order 1/N. Thus, if N 
is big, (AT)2/T^ is small and fluctuations may be expected to remain very small along 
the linear series up to or rather very near the stability limit at point D. In gravitational 
thermodynamics N is much smaller than in atomic physics; N ~ lO^'^ in galactic nuclei, 
10^~^ in globular clusters and 10^~^ in present day numerical simulations (Meylan and 
Heggie 1997). Fluctuations in these different cases are wildly different and some have 
measurable effects as we shall see. 

In this connection, it is worthwhile recalling Monaghan's (1978) application of the 
theory of hydrodynamic fluctuations to a self-gravitating gas. He showed that density 
fluctuations become large before the point of ordinary stability is reached. Statistical 
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mechanics confirms Monaghan's finding. 

In the present work we apply our own version of fluctuation theory in N body systems 
in statistical physics with gravitational forces (Parentani et al. 1995). The theory, first 
sketched in Okamoto et al. (1995), has been applied to evolutionary sequences of quasi- 
static equilibrium configurations of self-gravitating radiation in the presence of a black 
hole. Such systems exhibit a first order phase transition and fluctuation theory plays a 
useful role in explaining when superheated black holes or superheated radiation become 
unstable. 

Fluctuation theory plays a very different role in an evolutionary sequence of quasi- 
static isothermal spheres, for instance. These are in a state of local entropy maximum 
and the maximum is unique. The global entropy maximum is in fact infinite as shown by 
Antonov. The non-equilibrium entropy has thus a local maximum and a local minimum. It 
is qualitatively like the dashed curve in figure 2. Large enough fiuctuations in equilibrium 
configurations near marginal stability at point D may bring the system along the dashed 
line at a point where the entropy can grow towards the global maximum, at infinity. 
Our aim is to find a limit where global temperature fluctuations become large enough to 
put the system out of equilibrium. Such a limit exists if the mean square amplitude of 
temperature fluctuations grows exponentially. This is the case in general and as we shall 
see it is certainly so in isothermal spheres. 

Our theory of fluctuations applies not only to isothermal spheres but to any ther- 
modynamic equilibrium with long range forces which has a unique but local minimum of 
entropy and satisfy some additional conditions which we describe in detail below. 

In section 2 we review the theory of fluctuations that has been presented in Parentani 
et al. (1995) and in section 3 we apply the theory to isothermal spheres. As we shall 
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see, the stability limit induced by fluctuations depends exponentially on a power of N and 
the density contrast at which instability develops may be significantly lower than 709 for 
"small" N. The significance of these calculations is analyzed in section 4. 

2. Fluctuation theory with long range forces 

(i) Steepest descent expression for the entropy of a finite isolated system 

Consider a system of gravitationally interacting particles in a finite volume V, with 
total mass M and fixed energy E. The total Gibbs entropy is S = klnQ, where 0, is the 
volume of available phase space; it can be evaluated by a steepest descent method, the 
saddle point value giving the mean field entropy (HK77). Analysis of the quadratic fluctu- 
ations yields thermodynamical stability conditions and the next order term in calculating 
the entropy beyond the mean field approximation (see for instance Horwitz 1971). For a 
recent review on the subject of statistical mechanics of gravitating systems, see Padman- 
abhan (1990). In HK77's scheme Cl is given by a functional integral but one needs only 
consider continuous functions (Ginibre 1971) and, under general conditions (Courant and 
Hilbert 1953) which are usually met in physical problems, the continuous functions can be 
expanded in an eigenfunction series which converges absolutely and uniformly in the whole 
domain of existence. Under such circumstances, the functional integral can be replaced 
by discrete integrations in terms of an infinite but denumerable set of variables say x*. ft 
assumes thus a form like this 




(2.1) 



In the steepest descent evaluation, w is extremized. 



dw 



T = 0. 



(2.2) 



dx'^ 



6 



w is then expanded in powers of (x* — X*) = Ax^ where M, E) is a solution of (2.2). 

We shall assume that the x*'s are chosen in such a way that at the point of extremum 
indicated by a sub-index e, the matrix of second derivatives of w is diagonal; thus 



The minus sign is for convenience^. We can always denumerate the Aj's in such a way that 

Ai < A2 < A3 Expansion of w in powers of Aa;* may thus be written as follows (if A 

is properly chosen) 



The Aj's are known as Poincare's coefficients of stability following Poincare's (1885) seminal 
analysis of stability. In writing (2.4) we have nowhere assumed that we deal with pure 
or not pure gravitational interactions. (2.4) is therefore quite general and applies to non- 
extensive systems. 

Equation (2.4) cries for qualifying restrictions: 

(1) The summation like the integral is in general not convergent (see HK77 for details). 
What makes them converge is a short distance cutoff in the gravitational interaction. The 
effect of the cutoff is to modify the higher values in the spectrum Aj which play no role in 
stability and fluctuation considerations (see below). The effect of short distance cutoffs on 
stability was analyzed in Aronson and Hansen (1972). 

(2) w must admit second order derivatives with respect to x*. We shall assume that 
w is twice differentiable also with respect to V, M and E in the vicinity of the maximum 
otherwise thermodynamic coefficients would not exist. In what follows we shall concentrate 

J By assuming that variables may be chosen such that (2.3) holds, we bypass of course 
the main technical problem of stability theory which is to find the Aj's... 




(2.3) 




(2.4) 
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on the E dependence of equilibrium configurations and take V and M fixed. Therefore we 
write S{E), (3{E) etc... 

(3) We shall now regard the smallest A^'s to be different from each other. This assump- 
tion is more than we shall need and is made for convenience; it simplifies explanations. In 
the paper by Horwitz and Katz of (1978) it was shown that in stable isothermal spheres 
Ai > and that at density contrasts near but greater than 709, Ai < but A2 > 0. An- 
other case in which a few of the smallest Aj have been shown to be different is in models 
for stars of cold catalyzed matter; this includes cold neutron stars (see in Harrison et al. 
1965). At each turning point of the linear series —iJ,{N) where 11 is the chemical potential§ 
H = dE/dN there is a change of stability at each of the three turning points calculated and 
each corresponds to one and only one square eigenfrequency (a;^) of dynamical perturba- 
tions changing sign. Degenerate spectra of Aj's are in fact rare in astrophysical applications 
though not unknown. Degeneracy corresponds to bifurcations of equilibria; see for instance 
the many bifurcations in liquid ellipsoids linear series (Chandrasekhar 1969). They usually 
correspond to an excess of symmetry which can be lifted with small arbitrary perturbations 
of the potential energy [see Thompson (1979), for Maclaurin ellipsoids in particular see 
Katz (1979)]. There are a few exceptions of physical interest in which bifurcations cannot 
be removed (see Arnold 1986). When all we are interested in is the passage of stability to 
instability it is enough for Ai to be different from A2. Some interesting results about the 
onset of instability may also be obtained without assuming that the spectrum of A^'s is 
non degenerate (Sorkin 1981). 
(ii) The standard limit of stability 

With (2.4), we now define a "temperature function" T{E;x^) or rather the more 
§ The curve does not appear in Harrison et al.; it is represented schematically in Katz 

(1981) and is based on Harrison et al.'s tables. 
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convenient inverse temperature or /? = 1/kT which has the same units as £^ 

''=BE = -klE^^ ^-liE + O. . ;3 + A.— Ax- + O.. (2.5) 

i i 

The extremal value oi (3 is (3 — 1/ kT, T is the global temperature of equilibrium. For small 
values of Aa;*'s 

dX^ 

A^^P-Pc^^X,-^Ax\ (2.6) 

i 

These formal fluctuations are related to fluctuations of temperature in a small part of 
the system in the usual sense when the rest of the system can be treated as a reservoir 
[see Parentani et al. (1995) for a detailed discussion of this point]. That, however, is not 
normally the case in gravitating systems. 

Since w is twice differentiable, derivatives of P with respect to and E exist; the 
derivative of P defined in (2.5) with respect to a certain a;', keeping E constant and all 
other a;* = X* is readily derived from (2.6) 

op dX' (dp\ 

-— = Aj— — - = -— no summation on t\ (2.7) 

dx' ' dE \dx'J^ ^ ^ ^ ' 

and the derivative of $ with respect to E, keeping all = X* is obtained from (2.5) 

dE-dE-2^^\liE) ^^-^^ 

Now replace dX^/dE in (2.8) in terms of {d(3/dx'')e as given by (2.7); then dP/dE can be 
written 

dl^d£_^]_fdlV 

dE dE ^ \ ydx' j ' ^ ' ' 

We now turn our attention to the stability limit. The entropy is a maximum and the 
X*(£')'s represent a stable thermodynamic mean field configuration when all the Poincare 
coefficients are positive, i.e. when Ai > 0. Consider thus a linear series P{E) of stable 
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equilibrium configurations. Any change of stabifity or change in the degree of instabihty is 
characterized by a change of sign in one of the Aj's. Thus when a series of stable equilibria 
becomes unstable, Ai becomes small and changes sign. Near instability, the sum in (2.9) 
is dominated by the 1/Ai term assuming (ji)(3/dx^^ ^ 0. If (ji)(3/dx^^ = we can get 
rid of the coincidence of zeros by adding a small perturbation Sw{x'^) to w, for instance 
ex^E with e small enough. Then (^8(3/ dx^^ = e. This is the normal way bifurcations are 
removed in ordinary stability analysis (Hunt 1977)|. Small perturbations will generally get 
rid of such degeneracies (but not always). We assume that (^8(3/ dx^^ ^ or, if necessary, 
is non zero with a small perturbation of to . If the sum in (2.9) is indeed dominated by 
the first term which diverges when Ai = 0, and if both dP/dE and {dp/dx'^)e are finite, 
then for Ai 0, d(3/dE must diverge and it follows from (2.9) that 

, V 

dE Ai \ dx^ 

or that 



(2.10a) 



dE f dpV , 

This is Poincare's turning point property of E{P): the energy is maximum or minimum 
when a change of stability takes place. With (2.10) we can calculate the heat capacity Cy 
in stable configurations near instability: 

Thus Cy < for a stable equilibrium (Ai > 0) and Cy > for an unstable one. This 
property is true near a turning point. It is interesting to remind the reader that Cy < 



J Notice however that we do not know what the x*'s are in general to begin with (see 
footnote p. 7). For details on this point see Katz (1980). 
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in stable systems is a property of finite non-extensive microcanonical ensembles in thermal 
equilibrium and not of gravitating systems only. Notice that a microcanonical ensemble is 
here stable where a canonical ensembles is unstable {Cy < 0). 
(iii) Mean quadratic fluctuations of temperature and their probability 

Consider now a stable equilibrium configuration near a point of instability where 
Ai(> 0) — > and Aj > for z > 2. All terms of the quadratic sum in (2.4) except the 
one with Ai may be integrated out in (2.1); they are by definition strongly stable if they 
have strongly negative exponentials which is generally the case when N » 1. Thus, near 
instability, we are interested in the terms of w that are not (strongly) stable, i.e. in 

^ ^ _ iX^iAx'f (2.12) 

with a slightly renormalized factor A in (2.1). We may replace Ai in (2.12) by its expression 
given in (2.10b). Near instability, with all other 's integrated out, the (new) temperature 
function /? calculated from w{E;x^) reduces to 

P = P + Apc^P+(^^^ Ax\ (2.13) 
If we now replace Ax^ in (2.12) in terms of Ap given by (2.13) we obtain 



Equation (2.14) is of the same form as equation (110.3) of Landau and Lifshitz's (1985) 
chapter on (non-quantum) fiuctuation theory. We can now use the standard arguments 
of fiuctuation theory and say that the probability dW for a fiuctuation of P in the range, 
(3 + A(3 and (5 -\- A(3 -\- d(3, is proportional to exp(w — S); thus, the properly normalized 
dW is given by 

2 dp ^ 



dW = \ - — 7-5- exp 



27r dp 
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dp (2.15) 



From (2.15) it follows that the mean quadratic fluctuations of the global temperature, as 
a function of the equilibrium parameters, is given by 



The analogy with classical fluctuation theory is flagrant except for the sign of Cy- Notice 
that our analysis does not use small subsystems treating the rest of the system as a heat 
bath. Small subsystems have no well defined energies independently of the big system 
when there are long range forces. The intriguing difference of signs of Cy and comparison 
with familiar results is discussed in detail in Parentani et al. (1995). Before getting any 
further it is good to remember that (2.15) and (2.16) are valid only close to the point of 



marginal stability, point D, where dE/d(3 — 0, Cy — and A/?^ = oo. The derivative of 
E{j3) near point D may be approximated by the lowest order in a Taylor expansion: 

and with (2.17) the exponential factor in dW can be written 



exp 



(2.18) 



There appears here two dimensionless quantities which deserve a special symbol, say 

AU = ^—^ and = ^ (2.19) 

PD PD 

AU depends on equilibrium parameters only, while AV is a fluctuation in units of Pd- 
The term in (2.18) that contains second order derivatives of E is non-dimensional. The 
virial theorem for particles in a container tells us that the total kinetic energy and the total 
potential energy are of the same order of magnitude. Thus the total energy E is itself of 
the same order of magnitude than the total kinetic energy. On the other hand kT = 1/(3 
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is of the order of magnitude of the kinetic energy of one particle. Therefore the first factor 
in 2.18, which "goes hke (3E" must be of order N 

where the number Fd must be of order 1 « N. In terms of AC/, AV^ and Tr,, the 
probabihty dW defined in (2.15) can be written 



dW 



^^NAUexp -^roNAUiAVf dV. (2.21) 



The expression for the mean quadratic fluctuations (2.16) can be rewritten in terms of C/'s 
and y's in the foUowing form 



(AF)2 = {ToNAUy^ (2.22) 

The rate at which a small fluctuation of temperature disappears is, by virtue of the 
fluctuation-dissipation theorem, the inverse of the time it takes to return to equilibrium, 
1 / trei ■ This relaxation time is proportional to some power of N (see for instance Binney and 
Tremaine about stellar systems). But dW is a power of e~^^^ and if AU is not very small 
(say AU ~ 0.1), the exponent becomes extremely small for large A^( say N > 100). The 
probability per unit time dW/dt will thus also be dominated by the negative exponential 
of dW. 

(v) The destabilizing effect of fluctuations near a local maximum of entropy 

Let us now determine the role of global temperature fluctuations near marginal sta- 
bility. Consider a sequence of quasi-equilibrium conflgurations evolving with a relative 
rate of change in volume V{dt/dV) of a few trei- The quasi-static linear series has a local 
entropy maximum and a local entropy minimum (see flgure 2). The evolutionary sequence 
approaches the turning point D through a succession of quasi-equilibria. However, as point 
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D is approached, the mean quadratic fluctuations tend to infinity as can been seen from 
(2.16) since Cy and the system is unlikely to remain in thermal equilibrium. So how 
close to point D will the system survive in quasi-equilibrium? Let us try to define such a 
point, call it C, and find a posteriori if it has any physical sense. 

Let Tc be the temperature at point C and (AT')^^ its mean quadratic fluctuation. We 
define point C as one at which (see figure 2) the mean quadratic fluctuation is just equal 
to (Te — Tc)^. For smaller mean fluctuations big real fluctuations have only a negligible 
probability to put the system in a state of entropy at the left hand side of the minimum 
of entropy in figure 2 (in one relaxation period). If mean quadratic fiuctuations are bigger 
than {Te — T^)^ the probability for a real big fiuctuation in a relaxation period increases 
also. Big fiuctuations may bring the system in a state in which the entropy can increase 
indefinitely and if it can it will so that thermal equilibrium will be lost. Point C defines a 
limit between these different behaviors; at C 

^{^TYc ^Te-Tcc:^ 2{Td - Tc) (2.23) 

The last quasi-equality comes from a quadratic approximation near the horizontal tangent. 
In terms of t/'s and F's defined in (2.19), (2.23) can be written 

[W?]c = [mU)c?. (2.24) 

and if (2.24) holds then, following (2.22), 

{AU)c = i'iTDN)-^/^ (2.25) 

This equality defines AUc and thus Pc and also Tc- At point C the exponent in (2.18) is 
just equal to 

'\AUc) - '\Pd-Pc) ~ '\To-Tc) ■ ^'-''^ 
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(2.25) gives the following value for 13c 

Pc = Pd[1 + (AC/)c] = Pd[1 + i4TDN)-^/% (2.27a) 

The corresponding value of the energy Ec can be estimated by expanding E(P) near point 
D in powers of Pc — Pd, remembering that {dE/dl3)D = 0. Taking account of (2.20), we 
have, to second order in /3c — (3d, using also (2.27a): 



Ec^Ed 



1+ (IV2)^^-V3- 



4(/?E/iV)D 

{(3c, Ec) are the coordinates of point C. 



(2.276) 



3. Applications to Isothermal spheres 

Isothermal spheres have been studied in great detail by many people. We refer to Binney 
and Tremaine for a modern presentation of the theory of isothermal spheres. With a 
Maxwellian distribution of energy per particle, the density p(r) depends exponentially 
on the gravitational potential p ~ e~^^. This explains why Newton's equation can be 
written in the form [see equation (371) in Chandrasekhar (1934) or (4.15b) in Binney and 
Tremaine] 



dr \ dr 

The equation can be replaced by a set of two first-order diff'erential equations (Emden 
1907). Here we use the (it, v) variables introduced by Chandrasekhar; first define the mass 
within a radius r as //(r) 

11= [ ATTpr'^dr. (3.2) 

Then introduce the variables 

47r pv^ 2 

u = , V = 4:TTGmPpr and C = Inr. (3-3) 



1^ 
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Equation (3.1) is equivalent to the following pair of equations for u and v (Chan- 

drasekhar's equations (404) and (405)) 

du , . , 
— = u{3-u)-v, (3.4) 

— = 2v-- . 3.5 
dC u 

Tables of solutions of these equations have been published by Emden (1907) and Chan- 
drasekhar and Wares (1949). The total energy E of isothermal spheres is calculable in 
terms oi ub and vb-, the boundary values at r = i? of it and v. The non-dimensional total 
energy and inverse temperature in figure 1 are the following functions oi ub and vb', 

and from (3.3) alone, since //(i?) = hb = M = Nm, 

Reciprocally, in terms of E* and b 

UB = bE* + ^, VB = b(^bE* + ^y (3.8) 

From (3.8), (3.7) and from (3.4), (3.5) at r = i? we obtain a differential equation for 
the linear series E*{b) represented in figure 1: 

db b^ \ ^ bE* + \) ^ ^ 

Equation (3.9) can be integrated directly; however the parametric form (3.4)- (3.5) deduced 
from (3.6)-(3.7) is more convenient for numerical integration; with C,b = Ini? we have 

^ = -- [2{bE* + \f + (6 - l){bE* + i) + 6 - 2] , (3.10) 
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h{hE* + \). (3.11) 



The density contrast TZ is obtained by integrating twice equation (3.1). In terms of tt, v 
and C, a first integration gives 

^ = --• (3.12) 
Witfi (3.7), (3.11) and (3.12) it then foUows that 

^ = 4S = exp /' f ; , (3.13) 

Figure 1 has been obtained by integrating (3.10) and (3.11) setting C = at a point with 
ho = 1.549 and Eq ~ 0.327 calculated from Emden's tables. 

Let us now evaluate the quantities related to the stability limit triggered by fluctua- 
tions. Point D, has coordinates E^^ and a density contrast TZd all equal to 

ho = 2.03, E*ij = -0.335, Ud = 709. (3.14) 

The number Vd associated with point D, and defined in (2.20) can be derived from (3.9); 
it is non-dimensional and can be written using (3.6) and (3.7) like this 



Thus, following (2.27), the coordinates of the stability limit due to fluctuations are 

he ^hD{l + Abc/bD) = 2.03(1 + 0.297V"^/^) (3.16a) 
E^ ~ E*o(l + AE^/E*^) = -0.335(1 - 0.6SN-^/^) (3.166) 
at which point the density contrast, according to (3.13), is 

Thus, taking account of (3.16) we obtain 

7^c ^ 709 • exp (-3.30A^"^/^) (3.18) 
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4. Remarks on these results and some observational implications 

The Table gives the stabihty hmit induced by fluctuations for values of N in the range 
10 < A?" < oo. In addition to the coordinates {hc,E^) and the corresponding values of 
density contrasts TZci the table also provides relative changes in temperature Abe /bo and 
energy AE^/Efj which give indications about the validity of the linear approximation 
near point D. Notice that relative corrections of be and Eq do not exceed 15% for N as 
small as 10. The relative reduction of the density contrasts ATZc/TZd, on the other hand 
changes significantly. 

The numbers have a simple interpretation: in isolated sphere in slow quasi-static 
evolution towards higher and higher density contrast as described in the Introduction the 
gravitational catastrophe or core collapse appears at lower density contrasts than 709. The 
change from TZd to TZc is very small for N = lO^'"^ but quite significant for A?" ~ 10^ and 
becomes drastic at lower values of N. 

A gravothermal catastrophy, as is well known, is not like an avalanche; the central 
parts of the stellar system gets hotter while the outer parts are left behind [see Lynden-Bell 
(1999) for a recent revue on gravothermal catastrophy]. The instability induces a change 
in the evolution which accelerates progressively. Our point here is to note that the change 
does not happen at point D but rather earlier, at point C. 

How sharply is point C defined? To appreciate the sharpness of this new stability limit, 
consider the probability distribution dW at an inverse (non-dimensional) temperature 
b> be- Consider also a fluctuation (b — b) (refere to flgure 1 again) as big as 2(6 — ^d), big 
enough to put the system in a state of ever growing entropy. The probability distribution 
for such a fluctuation deflned by (2.21) can here be written in terms of 6's rather than t/'s; 
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we obtain for dW: 

dW = .465iV^/V^/2 ^-y3)dV (4.1) 



where 



Thus for a given r 



b = bc + {r-l){bc-bD). (4.3) 



For AT = 10^ (or N = 10^), at point b ~ 1.01 (1.12) higher than be with a factor r = 3, 
the exponential factor exp (— |r^) ~ 10~^ aheady. At 6 ~ 1.15 (1.18) higher than be with 
a factor r = 4 , exp ( — |r"^) ~ lO"-*^^. Thus, the probabihty (per unit time) for fluctuations 
to be big enough to render the system unstable changes eery steeply for small departures 
from be- In this respect, point C is rather sharply defined . 

The gravothermal catastrophe must thus take place at lower density contrasts because 
of the fluctuations and this should show up in A'"-body calculations with isothermal spheres 
evolving through quasi-static configurations. Possible implications for real globular clusters 
are as follows. Galactic globular clusters are of course not closed isothermal spheres. They 
are correctly modeled by Michie (1963)-King (1966) models which have truncated Gaussian 
distributions and are not in thermodynamic equilibrium. However to the extent that our 
results for isothermal spheres have any bearing on the behaviour of observed globular 
clusters, we speculate that a greater number of globular clusters appearing in the tables 
of Trager et al. (1993) than indicated there may actually be in an advanced stage of core 
collapse because most observed globular clusters whose parameter allows fitting them to 
energy truncated gaussian distributions are close to marginal stability, that is to point D 
(Katz 1980). 
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TABLE 



N 


be 


-E*c 


7^c 


Abe 
bo 




Eh 




A7^c 


CO 


2.03 


0.335 


709 












1 


10^ 


2.04 


0.335 


686 


.003 


6.0 


• 10- 


-5 


0.97 


10^ 


2.04 


0.335 


660 


0.006 


3.0 


• 10- 


-4 


0.93 


10^ 


2.06 


0.334 


608 


0.014 


1.4 


• 10- 


-3 


0.86 


10^ 


2.09 


0.333 


510 


0.029 


6.3 


• 10- 


-3 


0.72 


102 


2.16 


0.325 


348 


0.063 


3.0 


• 10- 


-3 


0.49 


10^ 


2.30 


0.289 


154 


0.140 


1.4 


■ 10- 


~1 


0.22 



Figure Captions 



Figure 1 

The figure represents the dimensionless expression of the total energy E* = 
RE / GN'^m? of isothermal spheres as a function of a dimensionless inverse 
temperature h = GNm? / RkT , where k is Boltzmann's constant, and G the 
gravitational constant. Numbers along the line represent density contrasts 
IZ = p{0)/p{R), p(0) is the density at the centre, p{R) is that at the boundary. 
TZ grows continuously in the direction of the arrow. 

Figure 2 



The full line represents the schematic behaviour of the linear series E{i3) 
near point D with corresponding points C and E of equal energy Ec- The 
dashed line is a schematic curve representing w{Ec; X ^ clS db function of x , 
> 1) =constant. 



